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Abstract 

A new analytic treatment of three-dimensional homogeneous Bose and Fermi gases in 
the unitary limit of negative infinite scattering length is presented, based on the S-matrix 
approach to statistical mechanics we recently developed. The unitary limit occurs at a 
fixed point of the renormalization group with dynamical exponent z = 2 where the S- 
matrix equals —1. For fermions we find Tc/Tp ~ 0.1. For bosons we present evidence that 
the gas does not collapse, but rather has a critical point that is a strongly interacting form 
of Bose-Einstein condensation. This bosonic critical point occurs at nXj. ~ 1.3 where n 
is the density and At the thermal wavelength, which is lower than the ideal gas value of 
C(3/2) = 2.61. 
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I. INTRODUCTION 



In the so-called unitary limit of a quantum Bose or Fermi gas, the scattering 
length a diverges. This occurs at a fixed point of the renormalization group, thus 
these systems provide interesting examples of interacting, scale-invariant theories 
with dynamical exponent z = 2, i.e. non-relativistic. They can be realized experi- 
mentally by tuning the scattering length to ±00 using a Feshbach resonance. (See for 
instance l2| and references therein.) They are also thought to occur at the sur 



ace 



of neutron stars. These systems have also attracted much theoretical interest js- 17 1. 
There have even been some proposals to use the AdS/CFT correspondence to learn 



about these models 



18|-l21|. 



Because of the scale-invariance, the only length scales in the problem are based 
on the density n^/*^ where d is the spatial dimension, and the thermal wavelength 
Ay = -^Itx jmT . Equivalently, the only energy scales are the chemical potential 
and the temperature T. The problem is challenging since there is no small paramater 
to expand in such as no?' . Any possible critical point must occur at a specific value 
of X = yu/T. This can be translated into universal values for ricAi^, or for fermions 
universal values for T^/Tp where ep = ksTp is the Fermi energy. For instance the 
critical point of an ideal Bose gas is the simplest example, where UcXp = C(3/2) = 
2.61. 



The present work is the sequel to 22|, where we used the S-matrix based formu- 
lation of the quantum statistical mechanics developed in {23, 24 1. This approach is 
very well-suited to the problem because in the unitary limit the S-matrix S = —1, 
and kernels in the integral equations simplify. In fact, this approach can be used to 
develop an expansion in 1/a. The main formulas for the 2 and 3 dimensional cases of 
both bosons and fermions were presented, however only the 2-dimensional case was 
analyzed in detail in 22|. Here we analyze the 3-dimensional case. 
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The models considered are the simplest models of non-relativistic bosons or 
fermions with quartic interactions. The bosonic model is defined by the action for a 
complex scalar field cf). 

s = J d'^dt i^<p^d,<p - ^ - j (1) 

For fermions, due to the fermionic statistics, one needs at least a 2-component field 

s = I d'^dt ( J2 - - (2) 

In both cases, positive g corresponds to repulsive interactions. The bosonic theory 
only has a f/(l) symmetry. The fermionic theory on the other hand has the much 



larger SO (5) symmetry. This is evident from the work 321] which considered a rel 



ativistic version, since the same arguments apply to a non-relativistic kinetic term. 



This is also clear from the work 17|] which considered an iV-component version with 
Sp(2N) symmetry, and noting that Sp(4) = S0(5). 

The interplay between the scattering length, the bound state, and the renor- 
malization group fixed point was discussed in detail from the point of view of the 
S-matrix in 22|. In 3 spatial dimensions the fixed point occurs at negative coupling 

= —An'^/mA, where A is an ultra-violet cut-off. For g less than g^,, there is a 
bound state that can Bose-Einstein condense (BEG), and for the fermionic case, this 
is referred to as the BEG side. As is approached from this side, the scattering 
length goes to +oo. The bound state disappears at g^. When g approaches g^ from 
above, the scattering length goes to — oo, and this is referred to as the BGS side. In 
this paper we work on the BGS side of the cross-over since the bound state does not 
have to be incorporated into the thermodynamics. On this side the interactions are 
effectively attractive. 

Theoretical studies have mainly focussed on the fermionic case, and for the most 



part at zero temperature, which is appropriate for a large Fermi energy. The bosonic 
case has been less studied, since a homogeneous bosonic gas with attractive interac- 
tions is thought to be unstable against mechanical collapse, and the collapse occurs 
before any kind of BEC. The situation is actually different for harmonically trapped 



gases, where BEC can occur [25|. However studies of the homogeneous bosonic case 



were based on a small, negative scattering length 26|-|29|. and it is not clear that the 



conclusions reached there can be extrapolated to the unitary limit. Since the den- 



sity of collapse is proportional to l/a[27|, extrapolation to infinite scattering length 
suggests that the gas collapses at zero density, which seems unphysical, since the gas 
could in principle be stabilized at finite temperature by thermal pressure. One can 
also point out that in the van der Waals gas, the collapse is stabilized by a finite size 
of the atoms, which renders the compressibility finite. In the unitary limit, there is 
nothing to play such a role. In the sequel we will present evidence that the unitary 
Bose gas undergoes BEC when nX^ ^ 1.3. This lower value is consistent with the 
attractive interactions. We also estimate the critical exponent describing how the 
compressibility diverges at the critical point. 

In the next section we define the scaling functions that determine the free energy 
and density, and derive expressions for the energy and entropy per particle, specific 
heat per particle, and compressibility. In section III the formulation of the unitary 



limit in 22[ is summarized. The two-component fermion case is analyzed in section 
IV. Evidence is presented for a critical point with Tc/Tp ^ 0.1, which is consistent 
with lattice Monte Carlo simulations. Bosons are analyzed in section V, where 
we present evidence for BEC in this strongly interacting gas. Motivated by the 



conjectured lower bound 30| for the ratio of the viscosity to entropy density rj/s > 
h/AirkB for relativistic systems, we study this ratio for both the fermionic and bosonic 
cases in section VI. Our results for fermions are consistent with experiments, with 
rj/s > 4.72 times the conjectured lower bound. For bosons, this ratio is minimized 



at the critical point where rj/s > 1.26 times the bound. 



II. SCALING FUNCTIONS IN THE UNITARY LIMIT. 

The scale invariance in the unitary limit implies some universal scaling forms j^. 
In this section we define various scaling functions with a meaningful normalization 
relative to free particles. 

First consider a single species of bosonic or fermionic particle with mass m at 
chemical potential /i and temperature T. The free energy density has the form 

^ = -C(5/2) (^— j Tc(/x/r) (3) 

where the scaling function c is only a function of x = fi/T. {( is Riemann's zeta 
function.) The combination ^JmT jlix = l/A^, where is the thermal wavelength. 
For a single free boson or fermion: 

hm Cboson = 1, Cfermion = 1 ^ (free partlcles). (4) 

m/t^o 2v2 

It is also convenient to define the scaling function q, which is a measure of the 
quantum degeneracy, in terms of the density as follows: 

nX^ = q (5) 

The two scaling functions c and q are of course related since n = —dT/dfi, which 
leads to 

q = C(5/2)c' (6) 

where c' is the derivative of c with respect to x. Henceforth g' will always denote the 
derivative of g with respect to x. The expressions for c and q for free theories will 
be implicit in the next section. 
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Also of interest are several energy per particle scaling functions. At a renormal- 
ization group fixed point, the energy density is related to the free energy in the same 
way as for free particles: 

where V is the volume. For a free fermion, in the zero temperature limit, the energy 
per particle E/N — )■ The Fermi energy is 

ej, = _ 37rVv^) (8) 



m 

The above definition can also be used for bosons. Since ei? = /i in the zero tempera- 
ture free fermionic gas, this leads us to define the scaling function ^: 

5^ ^ 505/2) /ey/^ ^ 

In the limit T — 0, i.e. x — )■ oo, ^ — ?■ 1 for a free fermion. 

A different energy per particle scaling function, ^, is meaningful as x — )> 0: 

A^3(2v^-l)C(5/2)~ 
NT 2(2v^-2)C(3/2) 

With the above normalization ^ = 1 for a free fermion in the limit x — > 0. In terms 

of the above scaling functions: 

- C(3/2)(2^-2)c 
(2^2-1) « 

The entropy density is s = —dJ-'/dT, and the entropy per particle takes the form 

' C(5/2)f5£/i^^ (12) 



n \ 1 

Next consider the specific heat per particle at constant volume and particle num- 
ber, i.e. constant density. One needs dx/dT at constant density. Using the fact that 
n oc T^^'^q, at constant density q oc T^^/^. This gives 



The specific heat per particle is then: 



N N\dTj^y A \ q q' 
The isothermal compressibility is defined as 



T 



where the pressure p = —T. Since n = N/V and is kept fixed, 

fdn^^\ 1 q' 1 fniTV'^ q' , , 



Finally the equation of state can be expressed parametrically as follows. Given n 
and T, one uses eq. ([5]) to find x as a function of n, T. The pressure can then be 
written as 

^ / C(5/2)c(x(.,T)) \ 
^ I g(a;(n,T)) ^ ^ ^ 

In order to compare with numerical simulations and experiments, it will be useful 

to plot various quantities as a function of q or T jTp: 



T / 4 \ 



T, 



2/3 



F 



:i8) 



III. TWO-BODY SCATTERING APPROXIMATION 



The main features of the two-body scattering approximation developed in [2J] are 
the following. Consider again first a single component gas. The filling fractions, or 
occupation numbers, are parameterized in terms of a pseudo-energy £(k): 



which determine the density: 



^^^^ = i^io^ ^19) 



_ . rf^k 1 



where s = 1,-1 corresponds to bosons, fermions respectively and /3 = 1/T. The con- 
sistent summation of 2-body scattering leads to an integral equation for the pseudo- 
energy £(k). It is convenient to define the quantity: 

y{k) = e-/^(^M— (21) 
where = k^/2m. Then y satisfies the integral equation 

^(k) = l + /3y (^G(k-k')jy-^ (22) 



The free energy density is then 



(2vr) 



3 



-s logfl — se 



2 eP^ 



(23) 



The kernel has the following structure: 



G = -^\og{l + iXM) (24) 

where Ai is the scattering amplitude and X represents the available phase space 
for two-body scattering. The argument of the log can be identified as the S-matrix 
function. In the unitary limit, 

^ = i = Hk3i^- 

and the S-matrix equals —1. The kernel becomes 

where the — sign corresponds to g being just below the fixed point (7^,, where the 
scattering length a — > -|-oo on the BEC side, whereas the -|- sign corresponds to 
a —7- —00 on the BCS side. As explained in the Introduction, we work on the BCS 
side. 
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The angular integrals in eq. (122|) are easily performed. Defining the dimensionless 
variable k = k^/2mT, the integral equation becomes 

r°° r 1 Z 

yU) = l + 4: dn' e(K-K')J^+Q(K'-K) — — (27) 

Jo I i - szy[K') 

where z = e^/^ is the fugacity and 9(k) is the standard step function equal to 1 for 
K > 0, zero otherwise. 

Finally comparing with the definitions in the last section the scaling function for 
the density and free energy are 

g(x) = A Trf^^^Ml^ (28) 
Vtt Jo - sy{K)z 

and 

^ = r (1 - szy{K)e-^) - ^4^^^^4) (29) 

V7rC(5/2) Jo V 2e''-szy{K)J 

The ideal, free gas limit corresponds to y = 1 where q = sLi3/2{sz) and c = 

sLi^/2{sz)/({5/2), where Li is the polylogarithm. The BEC critical point of the 

ideal gas occurs at /i = 0, i.e. q = C(3/2). 

Consider now two-component fermions with the action ([2]). Here the phase space 

factor X is doubled and since G oc 1/X, the kernels have an extra 1/2: 

Cfcrmi = 2^bose (30) 

Due to the SU(2) symmetry, the two-component fermion reduces to two identical 
copies of the above 1-component expressions, with the modification fl30|) . 



IV. ANALYSIS OF FERMIONS 



Recall that for 2 component fermions, the 4 is replaced by 2 in eq. f l27|l . and by the 
SU(2) symmetry, the system reduces to two identical copies of the one-component 
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equations of the last sections; below we present results for a single component, it 
being implicit that the free energy and density are doubled. 

The integral equation for y{n), eq. ( 1271) . can be solved numerically by iteration. 
One first substitutes yo = 1 on the right hand side and this gives the approximation 
yi for y. One then substitutes yi on the right hand side to generate y2, etc. For 
regions of z where there are no critical points, this procedure converges rapidly, and 
as little as 5 iterations are needed. For fermions, as one approaches zero temperature, 
i.e. X large and positive, more iterations are needed for convergence. The following 
results are based on 50 iterations. 

When z 1, y ^ 1, and the properties of the free ideal gas are recovered, since the 
gas is very dilute. There are solutions to eq. fl27|) for all — oo < x < oo. {x = ii/T). 
The scaling function c, and it's comparision with a free theory, are shown in Figure 
[1] as a function of x. The corrections to the free theory become appreciable when 
a; > -2. At a; = 0: 

c(0)= 0.880, ^= 0.884, (31) 

compared to the free gas values of c(0) = 0.646 and E, = ^-■ 

The scaling function q for the density is shown as function of x in Figure [2l Note 
that the density in the interacting case is always higher than for a free gas, due to 
the attractive interactions. At x = 0, g(0) = 1.18, whereas for a free gas q = 0.765. 
At low temperatures and high densities, fi/T ^ 1, the occupation numbers resemble 
that of a degenerate Fermi g 3iS, clS shown in Figure [3l 

Whereas c and q are nearly featureless, other quantities seem to indicate a critical 
point, or phase transition, at large density. For instance, the entropy per particle 
decreases with decreasing temperature up to x < Xc ~ 11.2, as shown in Figure |H 
Beyond this point the entropy per particle has the unphysical behavior of increasing 
with temperature. A further indication that the region x > Xc is unphysical is that 



10 



c for fermions 
c 




FIG. 1: c{x) and its equivalent for a free theory as a function oi x = fi/T. 

the specific heat per particle becomes negative, as shown in Figure O When x <^0, 
Cy/N approaches the classical value 3/2. This leads us to suggest a phase transition, 
at X = Xc, corresponding to the critical temperature Tc/Tp ~ 0.1. As we will show, 
our analysis of the viscosity to entropy-density ratio suggests a higher Tc/Tp. There 
have been numerous estimates of Tc/Tp based on various approximation schemes. 



mainly using Monte Carlo methods on the lattice lOMlSl]. quoting results for T^/Tp 



between 0.05 and 0.23. The work [ll| puts an upper bound T^/Tp < 0.14, and the 
most recent results of Burovski et. al. quote T^Tp = 0.152(7). Our result is thus 
consistent with previous work. The equation of state at this point follows from eq. 
(HZD: 

p = 4.95nr (32) 

The energy per particle, normalized to the Fermi energy ep, i.e. E/Nep = 3^/5, 
and the entropy per particle, are shown in Figures EE as a function of T/Tp, where 
ksTp = ep. At high temperatures it matches that of a free Fermi gas, in agreement 
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q for fermions 




= ^l/T 



FIG. 2: q{x) and its equivalent for a free theory as a function of x = /-f/T. 

13l . Il5| . Note that there is no sign of pair- 



with the Monte Carlo simulations in 
breaking at T* /Tp = 0.5 predicted in [10|, and this also agrees with the Monte Carlo 
simulations. However at low temperatures in the vicinity of T^, the agreement is 
not as good. This suggests our approximation is breaking down for very large z, i.e. 
the limit of zero temperature. The same conclusion is reached by examining fi/ep, 
displayed in Figured since the zero temperature epsilon expansion and Monte Carlo 
give iJ,/eF ~ 0.4 - O.sQ, Q- 



V. ANALYSIS OF BOSONS 



For bosons we again solved the integral equation (1271) by iteration, starting from 
y = 1. Since the occupation numbers decay quickly as a function of k, we introduced 
a cut-off K < 10. For x less than approximately —2, the gas behaves nearly classically. 

The main feature of the solution to the integral equation is that for x > Xc = 



12 



Occupation number for fermions 

/ 




FIG. 3: The occupation numbers as a function of k for a; = 5, 10, 15. 

— 1.2741, there is no solution that is smoothly connected to the classical limit x — )■ 
— oo. Numerically, when there is no solution the iterative procedure fails to converge. 
The free energy scaling function is plotted in Figure [91 Note that c < 1, where c = 1 
is the free field value. We thus take the physical region to be x < Xc- We find strong 
evidence that the gas undergoes BEC at x = Xc- In Figure [101 we plot £:(k = 0) as a 
function of x, and ones sees that it goes to zero at Xc- This implies the occupation 
number / diverges at k = at this critical point. One clearly sees this behavior in 
Figure [iTl 

The compressibility is shown in Figure [121 ^ind diverges at Xc, again consistent 
with BEC. We thus conclude that there is a critical point at Xc which a strongly 
interacting, scale invariant version of the ideal BEC. In terms of the density, the 
critical point is: 

n^\% = 1.325, (/i/T = Xc = -1.2741) (33) 
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entropy per particle for fermions 

s/n 




FIG. 4: Entropy per fermionic particle as a function of x. 

The negative value of the chemical potential is consistent with the effectively attrac- 
tive interactions. The above should be compared with the ideal BEC of the free 
theory, where Xc = and n^-X^ = C(3/2) = 2.61, which is higher by a factor of 2. At 
the critical point the equation of state is 

p = O.SlSnT (34) 

compared to p = 0.514nT for the free case. (0.514 = (^(5/2)/C(3/2)). 

A critical exponent u characterizing the diverging compressibility can be defined 

as 

K~(T-T,)^^ (35) 

A log-log plot of the compressibility verses T — shows an approximately straight 
line, and we obtain u ^ 0.69. This should be compared with BEC in an ideal gas, 
where u ~ 1.0. Clearly the unitary gas version of BEC is in a different universality 
class. 
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FIG. 5: Specific heat per particle as a function of x for fermions. 
Energy per particle scaling function for fermions 

E 
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1.0 



T/Tf 



FIG. 6: Energy per particle normalized to ep as a function of T/Tp- 

The energy per particle scaling function ^ at the critical point is ^(xc) = 0.281 
compared to 0.453 = (2^2 - 2)7(272 - 1) for the free case. The entropy per particle 
and specific heat per particle are plotted in Figures fT3l [14] as a function of T/Tc. At 
large temperatures, as expected Cy/N = 3/2, i.e. the classical value. It increases as 
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FIG. 7: Entopy per particle as a function of T jTp. 
Chemical potential for fermions 




FIG. 8: Chemical potential normalized to ep as a function of T jTp. 

T is lowered, however in contrast to the ideal g cLS CclSG, it then begins to decrease as 
T approaches Tc. 
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c 

0.35 r 



c for bosons 



0.30 



0.25 



0.20 



ideal 



-1.38 -1.36 -1.34 -1.32 -1.30 -1.28 



X = fi/T 



FIG. 9: The free-energy scaling function c as a function of /u/T compared to the ideal gas 



case. 



VI. ENTROPY TO VISCOSITY RATIO 



Consider first a single component gas. In kinetic theory, the shear viscosity can 
be expressed as 



(36) 



where v is the average speed and £free is the mean free path. The mean free path is 
^free = 1 / {V2na) whcrc a is the total cross-section. (The a/2 comes from the ratio 
of the mean speed to the mean relative speed 34].) In the unitary limit the S-matrix 
5* = — 1, which implies the scattering amplitude in eq. f l2B]) . This leads to 

m^\M\^ 167r 



(37) 



4:71 |k|2 

where |fc| is the momentum of one of the particles in the center of mass frame, i.e 
|ki — = 2|k|. This gives 



V 



m V 
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(38) 




FIG. 10: The pseudo-energy e at k = as a function of x = /-f/T. 

Since the equation (I7l) is the same relation between the pressure and energy of a 
free gas, and the pressure is due to the kinetic energy, this imphes 

\mv' = ElN=^-^T (39) 

Since the entropy density s = —dT/dT, one finally has 

v = ^(!Lf' 1 (40) 

s 8C(5/2) VcV 5c/2-xc' ^ ' 

For two-component fermions, the available phase space X is doubled. Also, spin up 
particles only scatter with spin down. This implies 77 is 8 times the above expression. 
Since the entropy density is doubled, this implies that 77/5 is 4 times the expression 
eq. dlO]). 

The ratio t]/ s for fermions as a function of T/Tp is shown in Figure [151 and is 
in good agreement both quantitatively and qualitatively with the experimental data 
summarized in 37|]. The lowest value occurs at x = 2.33, which corresponds to 
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FIG. 11: The occupation number /(k) for x = —1.275 and Xc = —1.2741. 
T/Tp = 0.28, and 

^ > 4.72-4- (41) 

s Air Kb 

The experimental data has a minimum that is about 6 times this bound. In the free 
fermion theory the minimum occurs at fi/T ^ 2.3, which gives rj/s > 1 .2h/ AirkB- 

For bosons, the ration rj/s is plotted in Figure [16] as a function of T jT^. One sees 
that it has a minimum at the critical point, where 

1 > 1.26-^ (42) 
s AnkB 

Thus the bosonic gas at the unitary critical point is a more perfect fluid than that 
of fermions. On the other hand, the ideal Bose gas at the critical point has a lower 
value: 

rj v/37rC(5/2) h 

= 20C(3/2)3/^ = '-''a^b ^''^ 

ideal 
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Compressibility for bosons 
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FIG. 12: The compressibility function of /x/T. 

Entropy per particle for bosons 
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FIG. 13: The entropy per particle as a function of T/Tc 
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Specific heat per particle for bosons 
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FIG. 14: The specific heat per particle as a function of T/Tc 
Viscosity/entropy for attractive fermions 
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FIG. 15: The viscosity to entropy-density ratio as a function of T/Tp for fermions. The 
horizontal line is l/47r. 

VII. CONCLUSIONS 

We presented a novel analytic treatment of unitary Bose and Fermi gases at finite 
temperature and chemical potential using a new formulation of statistical mechanics 



21 



Viscosity to entropy ratio for bosons 




FIG. 16: The viscosity to entropy-density ratio as a function of T/Tc for bosons. The 
horizontal line is I/Att. 

based on the exact, zero temperature, 2-body scattering. Our results appear to be 
consistent with lattice Monte Carlo methods. All of the thermodynamic functions, 
such as entropy per particle, energy per particle, specific heat, compressibility, and 
viscosity are readily calculated once one numerically solves the integral equation for 
the pseudo-energy. 

For fermions, our 2-body approximation is good if the temperatures are not too 
low. We estimated T^Tp ^ 0.1, where the critical point occurs at fi/T ^ 11.2. For 
bosons we presented evidence for a strongly interacting version of BEC at the critical 
point ~ 1.3, corresponding to fi/T = —1.27. 
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